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ABSTRACT 

One of the most well-known tests of General Relativity (GR) results from combining 
measurements of the anomalous precession of the orbit of Mercury with a determina- 
tion of the gravitational quadrupole moment of the Sun J 2 . The latter can be done by 
inference from an integral relation between J 2 and the solar internal rotation. New ob- 
servational data of high quality obtained from the Solar Heliospheric Satellite (SoHO) 
and from the Global Oscillations Network Group (GONG), allow the determination 
of the internal rotation velocity of the Sun as a function of radius and latitude with 
unprecedented spatial resolution and accuracy. As a consequence, a number of global 
properties of the Sun can also be determined with much higher accuracy, notably the 
gravitational quadrupole moment of the Sun. The anomalous precession of the orbit of 
Mercury is primarily due to GR effects but there are classical corrections the largest 
of which is that due to J 2 ■ It is shown here that the data are currently consistent with 
the predictions of GR. 
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1 INTRODUCTION 

For observations that are well resolved in space and in time 
the oscillations of the Sun can be decomposed into its pulsa- 
tion eigenmodes, which are products of functions of radius 
and of spherical harmonic functions. Each mode, and there- 
fore each measured oscillation frequency, is uniquely identi- 
fied by three numbers : the radial order n, and the degree I 
and the azimuthal order m of the spherical harmonic. The 
solar rotation produces oscillation frequencies that are split 
into multiplets. The relationship between the mode frequen- 
cies and the rotation is : 



2tt 



2m 



dx d cos 9 K n i m (x,6)Q(x, 9) (1) 



where x = r/i? Q is the fractional radius, Rq is the ra- 
dius of the Sun, and 9 the colatitude. The K n i m are the 
mode kernels for rotation. The Global Oscillations Network 
Group (GONG) produces values for the splittings through 
their data-reduction pipeline which are then available for 
inversion of the above integral relation. The Solar Helio- 
spheric Satellite (SoHO) SOI/MDI instrument pipe-line gen- 
erally produces Ritzwoller-Lavely a-coefHcients (Ritzwoller 
& Lavely, 1991), instead of individual splittings. The rela- 
tion between these a-coefficients and the rotation rate is a 
linear integral equation very similar to (1) although with dif- 
ferent kernels. Explicit expressions for the kernels for both 



cases can be found in e.g. Pijpers (1997). 

Using these data it is possible to determine the internal 
rotation rate of the Sun using inverse techniques. Results of 
such inversions can be found in e.g. Thompson et al. (1996) 
and Schou et al. (1997). Apart from the resolved rotation 
rate, there are some global properties of the Sun of astro- 
physical interest, which are related to the internal rotation 
rate through integral equations. One of these quantities is 
the total angular momentum H of the Sun which is related 
to the internal rotation rate through : 



l l 

H = J Ax J dcos9m{x,9) 

-1 

with the moment of inertia kernel T 
X = 2tyR%px a {1 - cos 2 9) 



(2) 
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where p is the density inside the Sun. Another is the total 
kinetic energy T in rotation which is given by : 



T = I <lr I dcos# ^m 2 {x,9) 



(4) 



Since the total angular momentum is related linearly to the 
rotation rate Q, it is possible to construct the kernel T di- 
rectly from a linear combination of the individual model 
kernels (or a-coemcient kernels) using for instance the tech- 
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nique of Subtractive Optimally Localized Averages (SOLA) 
(cf. Pijpers and Thompson, 1992, 1994) as was done using 
GONG data by Pijpers (1998). This avoids the circuitous 
route of first determining the resolved rotation rate and then 
re-integrating. The reasons for doing this are that it can 
have better properties from the point of view propagation 
of the measurement errors, as well as avoiding systematic 
errors introduced at each computational step and it is com- 
putationally much less expensive. It is not possible to do 
the same for the kinetic energy since this is quadratic in 
the rotation rate. To determine T, one has to either follow 
the route of determining the resolved rotation rate, taking 
the square and then re-integrating, or one must make use 
of the second order splittings which are generally much less 
well-determined and are affected by physical effects other 
than rotation such as internal magnetic fields (cf. Gough & 
Thompson, 1990). 

Another quantity of particular interest is the gravita- 
tional quadrupole moment J2 of the Sun, caused by its flat- 
tening due to the rotation. The gravitational quadrupole 
moment J2 of the Sun is that component of the gravita- 
tional field corresponding to the second Legendre polyno- 
mial as a function of co-latitude in an expansion of the 
gravitational field on Legendre polynomials. It is related to 
the solar oblateness A4, the ellipticity of the visible solar 
disk, as J2 = | A^. The gravitational quadrupole moment of 
the Sun modifies the precession of the orbits of the planets. 
Therefore in using for instance the precession of the orbit of 
Mercury for testing the prediction from GR it is necessary 
to know J 2 . Expressions for the integral relation between 
J2 and the internal rotation rate of the Sun have been de- 
rived for special cases of a rotation rate dependent on the 
radius only, or on simple parameterizations with respect to 
latitude (cf. Gough, 1981, 1982 ; Ulrich & Hawkins, 1981). 
More general expressions have been given by Dziembowski 
& Goode (1992) who expand the rotation rate by projection 
onto Legendre polynomials. However it can be shown that 
this is a somewhat cumbersome approach and quite simple 
expressions can be found even for a general distribution of 
Q(x,0). 

In section 2, the integral relation between the gravita- 
tional quadrupole moment of the Sun and a general internal 
rotation rate is given. In section 3 the results are given of 
performing the direct inversion for H and the values for T 
and J2 obtained by taking the square of the resolved ro- 
tation rate and re-integrating, using data from GONG and 
using data from SOI/MDI on board SoHO. Conclusions are 
presented in section 4. 



2 THE GRAVITATIONAL QUADRUPOLE 
MOMENT 

The general expressions relating the various moments of the 
gravitational potential of rotating stars to their rotation rate 
have been given by Goldreich & Schubert (1968) and by 
Lebovitz (1970). These lead to what is essentially Clairaut- 
Legendre equations for the moments. For convenience the 
steps will be briefly repeated here. Starting point are Pois- 
son's equation which relates the gravitational potential to 
the density distribution : 

V 2 = -4ttGp , 



and the equation of motion : 
pVcj> = Vp - pfi(r, 6) 2 r sin dOf 



(6) 



where <j> is the gravitational potential, G is the constant of 
gravity, p and p are the gas density and pressure respectively, 
and Q is the rotation rate which is a function of radius r and 
co-latitude 6. G^is a unit vector perpendicular to the rotation 
axis. Writing equation (6) out in components yields : 

9<t> dp n 2 W x2 

d<t> dp 2 , ,2 

in which u = cos 6. Following the treatment of Goldreich 
& Schubert (1968) and Lebovitz (1970) for slowly rotating 
stars all quantities are described in terms of perturbations 
of the spherically symmetric non-rotating star, i.e. Q 2 is 
treated as a quantity of first order in a small parameter 
expansion. Subscripts refer to the non-rotating configura- 
tion, and 1 to the perturbed quantities. Collecting the first 
order terms in the perturbation analysis of equation (7) : 



d<)>i d(f>o _ dpi 
dr Pl dr dr 
d(pi _ dp 
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(8) 



Of interest for the quadrupole moment of the gravitational 
potential is the projection onto the Legendre polynomial 
Pi{u) = (3u 2 — l)/2. In the first equation of (8) all terms 
are multiplied by §P 2 («) and then integrated over u. The 
second equation would yield = since all its terms are 
odd in u. Therefore this equation is first integrated in u and 
then projected. In the following the subscripts 12 refer to the 
part of the first order perturbed quantities corresponding to 
these P2 Legendre polynomial projections. 



d<j>ii d</) dpi2 f 5 

p0 ~cV~ + pl2 ~dV = ~~d~F~ " p0 / 3 [ 2( ' 



>\i—pi2=por I du -P 2 (u) / dvvfl(r,v) 



(9) 



The double integral in the second equation can be re- written, 
using partial integration : 



dtt —P2(u) I dv vQ(r,v) = 



(5) 




1 (10) 
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The second equality of equation (9) can be used to eliminate 
P12 from the first of (9). After some rearranging the result 
is : 



P12 



dr 



i 



Por J du^(l~P 2 (u))P 2 (u)n(r,u) 2 ( n ) 



P12 = Po0i2 - por G(ty 
where Q is defined by : 
l 

G{Q) = ^Jdu(u- u 4 ) Q(r,uf 



(12) 



The relevant equations from the perturbed Poisson's equa- 
tion (5) are : 
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r 2 dr 



d 2 cj>i2 2 <90i 2 6 



(13) 



dr 2 r dr 



-012 = -4-KGpl 



in which p\2 can now be substituted using the first of equa- 
tions (11) : 

<9 2 0ia | 2 d<f>i2 _ 6_ = 47rr 2 LgPo 
dr 2 r dr r 2 M r 1 dr 
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dr 



[p r 2 g(Sl)} 



(14) 



-Por I du-{l-P2(u))P2(u)n(r,uY 



in which use has been made of the mass within radius r : 



r 

Mr = J dr' 
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Now define the linear differential operator £ 
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dr \ dr J 
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Mr ~dr 
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and a function / : 



f(r) 
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so that 0i2 is the solution of C(f)\ 2 = f(r). This equation can 
be solved used Green's functions. For r > Rq the density 
po = and therefore f(r) = 0. An exact solution is then 



0i2 = r .If another solution of Crp = is constructed 
which is regular at r = the general solution is : 



oo 



oi2(R) = J dz G(R,z)f(z) 



with the Green's function 

-3 



(18) 



G(R,z) 
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z 2 W{z 
j>{z)R- 
z 2 W(z 



0<R<z 



< z < R 



(19) 



where W(z) is the Wronskian of the solutions r 3 and tp : 

-3 

. =_*-« 

dz 



W(z) = 



z 
0' -3z~ 4 



= - z ^4-(z 3 ^) 



(20) 



Since r~ 3 is not a solution of Cip = for r < Rq this 
equation is only valid for R > Rq. Of interest here is the 
solution 0i2 at R = Rq. If in (18) R is replaced with Rq it 
is allowed to replace z 2 W(z) with RqW(Rq), which can be 
verified by substitution of (18) into (14). Since f(r) = for 
r > Rq, the expression for 0i2 (Rq) simplifies to : 



M^e) = -Re 



r — (r 0) 



dz0(z)/(z) (21) 



The solar oblateness is related to 0i2 as 

Substituting the expressions for 0i2 and f(r] 
9^Rl r A „ i - 1 

A 



(22) 
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in which the second equality is obtained by partial integra- 
tion, and the third is a re-arranging of terms making use of 
the definition (12) of Q. 

For an which is a function of r only, integration over u 
of the second term between square brackets in (23) is identi- 
cal to and the first term reduces to po-§^ jti^ ^j ^( r i u ) 2 ■ 
Equation (23) then reduces to equation (12) of Gough 
(1981). Once the density po(r) of the Sun is known, it is 
trivial to calculate the two-dimensional kernel : 



Normalized kernel J 
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(24) 
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Determining A^> is thus reduced to evaluating the two- 
dimensional integral : 



dr / du F(r, u)Q(r, uf 



(25) 



Direct inversion would have to make use of the second order 
splittings in an inverse problem, so the same route is followed 
as with the kinetic energy T : the resolved Q 2 is multiplied 
with T and integrated. 

Using the standard solar model S of Christensen- 
Dalsgaard (cf. Christensen-Dalsgaard et al., 1996) the ker- 
nels X and T were calculated, and normalized to have unit 
integral over the solar volume. Contour plots in one quadrant 
are shown in figure 1. The other quadrants can be obtained 
by reflection in the coordinate axes. 
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Figure 1. The normalized kernel I for determination of the total 
solar angular momentum and kinetic energy (top panel), and the 
kernel T for the determination of the gravitational quadrupolc 
moment of the Sun. Contour intervals are 0.1 for both kernels. 



3 RESULTS AND CONCLUSIONS 

Two independent data sets have been used to determine 
the total solar angular momentum, the total kinetic energy 
and the gravitational quadrupole moment. One dataset is in 
the form of splittings obtained with the earth-based GONG 
network of telescopes : 33169 splittings distributed over 542 
complete multiplets with 7 < I < 150 and 1.5 mHz < v < 
3.5 mHz gathered from GONG months 4 to 10. The other 
dataset is in the form of a-coefficients gathered from 144 
d out of the first 6 months of operation of the SOI/MDI 
instrument on board the SoHO satellite. The data consists 
of 414 multiplets with 1 < I < 250 and 1.0 mHz < v < 
4.2 mHz, and the odd a-coefficients up to at most (135 are 
available. 

The GONG data leads to the values : 
H d = [186.3 ± 2.4] 10 39 kg m 2 s _1 
Hi = [186.3 ± 3.7] 10 39 kg m 2 s" 1 
T = [245.5 ± 9.8] 10 33 kg m 2 s~ 2 
J 2 = [2.14 ±0.09] 10~ 7 
The MDI data leads to the values : 
H d = [192.3 ± 1.9] 10 39 kg m 2 s _1 
Hi = [192.9 ± 3.9] 10 39 kg m 2 s" 1 
T = [262.5 ± 10.] 10 33 kg m 2 s~ 2 
J 2 = [2.23 ± 0.09] 10~ 7 



(26) 



(27) 



The subscript d refers to a determination directly from the 
data using the freedom of the SOLA method to construct the 
kernel directly, the subscript i refers to the indirect method, 
which is re-integrating the resolved Q. T and J2 have been 
determined by re-integration only. Since the direct method 
should suffer much less from systematic effects, the value 
Hi is shown merely to demonstrate consistency between the 
two methods. Error weighted means for H d ,T, and J2 are : 

H = [190.0 ± 1.5] 10 39 kg m 2 s" 1 

T = [253.4 ± 7.2] 10 33 kg m 2 s~ 2 (28) 
J 2 = [2.18 ±0.06] 10~ 7 

This determination of J2 is entirely consistent with that of 
Paterno et al. (1996) who used direct oblateness measure- 
ments of the solar disk to infer the quadrupole moment. 

One of the most well-known tests of GR results from 
combining measurements of the precession of the orbit 
of Mercury (cf. Shapiro et al., 1976 ; Anderson et al., 
1987, 1991, 1992) with a determination of the gravitational 
quadrupole moment of the Sun J2. In the fully conservative 
parameterized post-newtonian (PPN) formalism, the pre- 
dicted advance A^o per orbital period of a planetary orbit 
with semi-major axis a and eccentricity e, after correcting 
for perturbations due to other planets, is : 



A0o 



6ttGM\ p 
a(l - e 2 )c 2 



(29) 
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where 



A 



-(2-/3 + 2 7 ) + 



-Ja 



(30) 



3 V ' " ' 2GMa(l-e 2 ) 

Here M and 7? are the mass and radius of the Sun, G is 
the gravitational constant, and c is the speed of light. The 
parameters /3 and 7 are the Eddington-Robertson parame- 
ters of the PPN formalism (cf. Misner et al., 1973), which 
in general relativity are equal to 1. For Mercury the above 
relation (29) reduces to : 

A(f> = 42.9794 A p "/century (31) 
and (30) is : 

A p = i(2-/3 + 2 7 )+2.96 10 3 x J 2 (32) 

Shapiro et al. (1976), using planetary radar ranging, found 
an anomalous precession for Mercury's orbit of 43.11 ± 0.21. 
Using radar and spacecraft ranging Anderson et al. (1987) 
found 42.92 ± 0.20 and an update (Anderson et al., 1991) 
gives the value 42. 94 ±0.20. The most recent result reported 
by Anderson et al. (1992) is 43.13±0.14. Combining the most 
recent value for the anomalous precession of Mercury's orbit 
and the value for J2 given above in equations (28) yields : 



-(2 - /3 + 27) = 1.003 ± 0.003 



(33) 



The error quoted here is entirely due to that in the planetary 
ranging data, since the error due to the uncertainty in J2 is 
two orders of magnitude smaller. 

In this paper it is thus demonstrated that the total solar 
angular momentum, its total kinetic energy in rotation, and 
the solar gravitational quadrupole moment can be deter- 
mined through inverting integral equations that are linear 
in the rotation rate Q or in its square, with known inte- 
gration kernels. The value of the gravitational quadrupole 
moment (28) when combined with planetary ranging data 
for the precession of the orbit of Mercury yields a value for 
the combined PPN formalism parameters (33) which is con- 
sistent with GR in which this combination is predicted to 
be exactly equal to unity. More stringent tests of GR using 
the orbit of Mercury rely on measuring its orbital precession 
with much greater precision. 



Recent Developments in Solar-System Tests of General 
Relativity, Proceedings of The Sixth Marcel Grossmann 
Meeting, Eds. H. Sato and T. Nakamura, World Scien- 
tific Publishing Co., Singapore 

Anderson J.D., Colombo G., Espsitio P.B., Lau E.L., Trager 
G.B., 1987, Icarus 71, 337 

Anderson J.D., Slade M.A., Jurgens R.F., Lau E.L., Newhall 
X X, Standish Jr. E.M., 1991, Pub. Astron. Soc. Aus- 
tralia 9, 324 

Christensen-Dalsgaard J., et al., 1996, Science 272, 1286 

Dziembowski W.A., Goode PR., 1992, ApJ 394, 670 

Goldreich P., Schubert G., 1968, ApJ 154, 1005 

Gough D.O., 1981, MNRAS 196, 731 

Gough D.O., 1982, Nature 298, 334 

Gough D.O, Thompson M., 1990, MNRAS 242, 25 

Lebovitz N.R., 1970, ApJ 160, 701 

Misner C, Thorne K.S., Wheeler J.A., 1973, 'Gravitation', 

Freeman, San Francisco, 1072, 1116 
Paterno L., Sofia S., DiMauro M.P., 1996, A&A 314, 940 
Pijpers F.P., 1997, A&A 326, 1235 

Pijpers F.P., 1998, Proc. IAU symposium 181, 'Sound- 
ing Solar and Stellar Interiors' Eds. J. Provost, F.-X. 
Schmider, Kluwer, Dordrecht, in press 
Pijpers F.P., Thompson M.J., 1992, A&A 262, L33 
Pijpers F.P., Thompson M.J., 1994, A&A 281, 231 
Pijpers F.P., Thompson M.J., 1996, MNRAS 279, 498 
Ritzwoller M.H., Lavely E.M., 1991, ApJ 369, 557 
Schou J., et al., 1997, ApJ submitted 

Shapiro I.I., Counselman C.C. Ill, King R.W., 1976, Phys. 

Rev. Lett., 36, 555 
Thompson M., et al., 1996, Science 272, 1300 
Ulrich R.K, Hawkins G.W., 1981, ApJ, 246, 985 



This paper has been produced using the Royal Astronomical 
Society/Blackwell Science TpX macros. 



ACKNOWLEDGMENTS 

The Theoretical Astrophysics Center is a collaboration be- 
tween Copenhagen University and Aarhus University and 
is funded by Danmarks Grundforskningsfonden. GONG is 
managed by NSO, a division of NOAO that is operated by 
the Association of Universities for Research in Astronomy 
under co-operative agreement with NSF. The GONG data 
were acquired by instruments operated by the BBSO, HAO, 
Learmonth, Udaipur, IAC, and CTIO The MDI project op- 
erating the SOI/MDI experiment on board the SoHO space- 
craft is supported by NASA contract NAG5-3077 at Stan- 
ford University. 



REFERENCES 



Anderson J.D., Campbell J.K., Jurgens R.F., Lau E.L., 
Newhall X X, Slade III M.A., Standish Jr. E.M., 1992, 



